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INTRODUCTION 


A line plotting routine has been developed using Langley's DI-3000 
graphics libraries. The program is written in Fortran 5 primarily using 
Graf maker subroutines (Ref. 1) and is general enough for practical use. 

Proqam Capabilities - The program allows multiple lines on a frame from 
multiple files with separate independent variable sets. Also, multiple 
frames are allowed within each run. Different linestyles and symbols are 
incorporated to distinguish data sets as well as full legend capability 
and titles. The option exists for automatic or manual scaling; linear, 
log or semi-log plots. The number of data points per line and the 
minimum and maximum values of the data need not be known to plot 
information. DI-3000 graphics allows greek and other fonts as well as 
subscripting and superscripting in the plot text (Ref. 2) . Combined with 
the Metafile translator (Ref. 3), graphic output has numerous formats. 

Program Limitations - Data to be plotted must be formatted such that 
within the data file, the independent data is located in the first 
column and dependent data is located in columns 2-6. Unless the data 
is to be read from the first two columns only (i.e. only one set of 
dependent data corresponds to a set of independent data) , the data 
must be in 20 character width columns. Log plots are limited to one 
cycle, and the minima and maxima must be specified. Currently, the 
program is limited to 7 lines per frame and 500 points per line. 


PLOT CONTROL FILE 

The program was written such that information required by the 
plotting routine is obtained from a "control file". The control file has 
the advantage of storing this information so that necessary changes can 
be made without repeating input of unchanged information. A control file 
has the following format: 


Line # 

Variables 

Description 

1 

NFRAME 

- Number of Frames 

2 

TITLE 1 

- First Line of Frame Title, Text 
Surrounded by $ signs (Delimeters) 

3 

TITLE2 

- Second Line of Frame Title, Text 
Surrounded by $'s (Delimeters) 

4 

LABELX 

- Label of X-Axis, Text Surrounded by $'s 

5 

LABELY 

- Label of Y-Axis, Text Surrounded by $ ' s 

6 

LEGEND 

- Legend Title, Text Surrounded by $'s 
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7 HFORM,HFW - Fortran Format Specification for X-axis Tic 

Labels, Width Value of Format Specification, 
Surround Format Specification with $'s, 

Line is Read in A10,I4 format Ex.- $(F3.1)$ 3 
Note: HFW is the same number as width value 

of Fortran Format Specification 

8 VFORM , VFW - Fortran Format Specification for Y-axis Tic 

Labels, Width Value of Format Specification, 
Surround Format Specification with $'s. 

Line is read in A10,I4 format Ex.- $(F4.2)$ 4 

Note: VFW is the same number as width value 

of Fortran Format Specification 


9 NLINES , HTICINC , VTICINC , I GRID, IHLOG , IVLOG , IINP 

NLINES - Number of Lines on Frame 

HTICINC - Increment Between Horizontal Tic Marks 

(HTICINC=0 , IINP=0 , and HMIN, HMAX must 
be specified for default tic marks) 

VTICINC - Increment Between Vertical Tic Marks 

(VTICINC=0 , IINP=0 , and VMIN, VMAX must 
be specified for default tic marks) 

I GRID = 0 - Grid at Major Tic Divisions 

= 1 - No Grid 

IHLOG, IVLOG - Horizontal, Vertical Axis Type 

= 1 - Linear Axis 
= 2 - Log Base 10 Axis 
= 3 - Log Base e Axis 

IINP = 0 - Program Will Expect Maxima and Minima on 

Next Input Line (IINP = 0 for Log Plots) 

= 1 - Program Generates Maxima and Minima for 
the Frame from Input Data 

[10] HMIN, HMAX, VMIN, VMAX 

- Maxima and Minima of Data. Input if IINP = 0 

10 [11] DFIL - Name of Data File Containing Data for Line, 

Independent and Dependent 

11 [12] NCOL, LANS, MARK, LAS 

NCOL = 2-6 - Column Number in Data File (DFIL) Containing 

Dependent Variable for Line. Data for Indep- 
endent Variable is Assumed to be in Column 1 

LANS = 0 - A Legend Entry for this Line Follows 

= 1 - No Legend Entry for this Line 
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MARK 


- Identification for Symbols to be Drawn 
= 0 - No Symbols, line only 
= 1 - Period " . " 

= 2 - "+" 

= 3 - "*" 

= 4 - "O" 

= 5 - "X" 

Note: If both MARK and LAS = 0 then a solid 

line will connect data points 

LAS = 0 - Draw Symbols Only 

= 1-7 - Linestyles (1-Solid, 2-Short Dash, 3-Long 

Dash, etc. Ref. 4 for Linestyle Information) 

[12] LEGN - Legend Text for this Line Surrounded by $'s 

Repeat lines 10 [11] through [12] NLINES times, then repeat lines 2 
through [12] NFRAME times. 


The data in the control file is read in free format except for the plot 
text and the format information (HFORM, VFORM) which must be enclosed by 
$ delimiters. The maximum length of the text strings is 50 characters 
for the titles, 30 for the legend title, 30 for each legend entry and 30 
for the axis labels. Filenames may not exceed 7 characters and need not 
be enclosed by delimiters. 


INTERACTIVE PROCEDURE 


For an interactive run of the line plotting routine, type: 


ATTACH , GMKCLIB, DI3CLIB/UN=LIBRARY . 

GET , LPLTRB/UN=8 2 6 0 3 5N . 

GET , CONTROL, DATAFI LE 1 , DATAFILE2 , DATAFILE3 
REWIND,*. 

LDSET , LI B=GMKCLI B/ DI 3 CLIB , MAP=N . 

LDSET , SUBST=DDNCAP-DD4014/MFNCAP-MFNODE/SSNCAP— SSDUMMY . 
LPLTRB . 

Where : 


CONTROL 
DATAFI LEI 

DD4014 


DDNCAP 

DI3CLIB 

GMKCLIB 


Plot Control File 

First Data File containing data to be plotted as 
specified in the control file 

Selected Device Driver (See section 5 of the on-line 
document "DINTRO" under UN=PVINFO for other device 
drivers - use XEDIT command L/5. DEVICE DRIVER/ 2 to 
find the beginning of this section. DD4014 is the 
Textronix 4014 device driver. 

Device Driver Capsule Library 
DI3000 Capsule Library 
Graphmaker Capsule Library 
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MFNCAP 
MFNODE 
LPLTRB 
SS DUMMY 
SSNCAP 


Metafile Node Capsule Library 

Metafile Driver 

Plotter Object Code 

Segment Storage Dummy Node 

Segment Storage Node Capsule Library 


The program then responds with: 

PLOT TO SCREEN ? (Y/N) 

- enter Y if you wish to view the plot on the terminal screen 

- enter N if you wish to create a Metafile only 

DO YOU WISH TO CREATE A METAFILE ? (Y/N) 

- enter Y if you wish to store the plot in the Metafile "DIMETA" 

- enter N if you wish to plot to the screen only 

ENTER CONTROL FILENAME: 

- enter the filename of the plot control file 

Plotting will then commence. If you created a Metafile, you may view the 
plot(s) again. To run the Metafile translator, type: 

ATTACH , DI 3 0 0 0 , MFDUMMY , SS DUMMY , DD4 014, DIMFT/UN=LI BRARY . 

REWIND, *. 

LDSET , LIB=DI3 000 , MAP=N . 

LOAD , MFDUMMY , DD4 0 14 , SS DUMMY . 

DIMFT. 

Where: 

DD4014 - Selected Device Driver 

DI3000 - DI3000 Graphics Library 

DIMFT - Metafile Translator Obect Code 
MFDUMMY - Metafile Dummy Node 
SSDUMMY - Segment Storage Dummy Node 

The Metafile translator responds with the prompt: 

MF> ? 

See Reference 3 for Metafile Commands. A sample run might look like: 

MF> ? S MF 1 DIMETA - sets the metafile to be viewed 

MF> ? D P 1 MF 1 - draws the first picture of Metafile 1 on the 

entire screen. 

Each time the plotting routine is run, "DIMETA" is overwritten, so 
the Metafile should be renamed and saved. To rerun the program again in 
the same session, only the last four commands of the procedure need be 
repeated. It is suggested that these commands be executed from a 
Procedure File (Ref. 4, Section 4) . To do this, simply type the 
procedure commands in XEDIT mode with the first line containing the 
statement: ". PROC, filename" 
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For a hard copy of the plot(s) , a Metafile must be created, and 
the following form of a Batch file submitted: 

/JOB 

jobname,T200. 

USER, . . . 

CHARGE , . . . 

DELIVER 

GET, metafile. 

PLOT. device* 

CONT. //Instructions to 
CONT . Operator/ / 

DAYFILE, jobok. 

REPLACE , j obok . 

EXIT. 

DAYFILE, badj ob 
REPLACE, bad job 
/EOR 

S MF 1 metafile |_ Metafile Translator Commands (Ref. 3) 

D P 1 MF 1 | 

/EOF 


♦Device = Plot Processors: 
ex. PLOT. CAL, 11 - Calcomp 11" Plotter 

PLOT. VAR, R - Varian Roll Plotter 

PLOT. VAR, F - Varian Fanfold Plotter 

To submit this batch file, type: 


SUBMIT, filename, B 


BATCH PROCEDURE 

Submit File for Batch Run of Line Plotter - For a large number of plots, 
it is advantageous to run the program in batch mode rather than 
interactive. To accomplish this, the following form of a batch file must 
be submitted: 


/JOB 

SUBMIT, T2 00. 

USER, .... 

CHARGE, . . . , LRC. 

DELIVER 

GET, control, datafilel,datafile2,datafile3, . . . 

ATTACH , DI 3 CLI B , GMKCLI B/UN=LI BRARY . 

GET , LPLTRB/UN=8 2 6 0 3 5N . 

LDSET , LI B=DI 3 CLI B/GMKCLI B , MAP=N . 

LDSET , SUBST=DDNCAP-DD4 014 /MFNCAP-MFNODE/SSNCAP-SSDUMMY . 

LPLTRB . 

REPLACE, DIMETA=metafile. - rename metafile 

PLOT. device - plot processors - Varian, Calcomp 
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DAYFILE, jobok. 

REPLACE , j obok . 

EXIT. 

DAYFILE, bad j ob. 

REPLACE , bad j ob . 

/EOR 

N 

Y 

control 

/EOR 

S MF 1 metafile |_ Metafile 

D P 1 MF 1 | 

/EOF 


To submit this batch file, type : 
SUBMIT, filename, B 


- leave one space before these 

- plot control file 
Translator Commands (Ref. 3) 
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Examples - Several examples are provided to show program capability. 
In each example, the plot control file is shown immediately preceeding 
the corresponding plot(s). 


Example 1 

This example shows the format of a typical control file: 


1 

$$ 

$$ 

$X/C$ 

$-C[BSUB] [BLC] P[ESUB] [ELC]$ 
$$ 

$ (F3 . 1) $ 3 

$ (F3 . 1) $ 3 

2. . 2 . .5. 1.1. 1.0 
—.2,1. ,—.5,1.5 
CP75 

2 . 0 . 0. 1 

$WITHOUT WAVE$ 

CP63 
2 , 0 , 0 , 2 
$WITH WAVE$ 


Number of Frames 
First Title, blank 
Second Title, blank 
X - Axis Label 
Y - Axis Label 
Legend Title, Blank 

Horizontal Label Format Specification 

Vertical Label Format Specif ication 

Line, Tics, Grid and Axis Information 

Maxima and Minima 

Data File for First Line 

Column, Legend and Line Information 

Legend Entry For First Line 

Data File for Second Line 

Column, Legend and Line Information 

Legend Entry For Second Line 


Note: Text Commands in Brackets cause subcripting, superscripting and 

font changes. For a description of embedded text commands, see Reference 
2, pp. 102-104. For a list of available fonts, see Reference 2, pp. 
93-96. 
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Example 2 


This example illustrates the use of embedded text commands (i.e. 
subscripting) and automatic scaling. Also, the example shows the first 
two linestyles. 


1 

$$ 

$$ 

$x/c$ 

$-C[BSUB] [BLC]P[ESUB] [ELC] 

$$ 

$ (F3 . 1) $ 3 

$ (F3 . 1) $ 3 

2, .2, .5, 1,1,1, 1 
cfilel 
2 , 1 , 0,1 
cf ile2 
2 , 1 , 0, 2 


- a lower case "p" will be subscripted 


- automatically find maxima and minima 
from the following data sets 

- no legend information 
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Example 3 


This example shows plot text capabilities such as first and second 
titles, and legend information. Also, manual scaling is shown. 


1 

$EFFECT OF MACH NUMBER$ 
$ONE-FOURTH SINE STEP$ 

$X/C$ 

$-C[BLC] [ BSUB ] P [ ESUB ] [ELC]$ 

$$ 

$ (F3 . 1) $ 3 

$ (F3 . 1) $ 3 

3, .2, .5, 1,1, 1,0 
—.2,1. , — .5,2.0 
RAMPCP1 
2, 0,0,1 
$M = 0.55$ 

RAMPCP2 
2 , 0 , 0 , 2 
$M = 0.65$ 

RAMPCP3 
2 , 0 , 0, 3 
$M = 0.75$ 


- two title lines 


- a lower case "p" will be subscripted 


- three lines, .2 between each horizon- 
tal tic .5 between each vertical 
tic, no grid, linear axes, and the 
minima and maxima follow. 
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Example 4 


This example shows the use of symbols and a grid. Also, all data for 
this plot comes from one data file. 


1 

$$ 

$M '= 0.750 ALT = 38000 FT.$ 
$X/C$ 

$-C[BSUB]P[ESUB]$ 

$$ 

$ (F4 . 2 ) $ 4 

$ (F5 . 2 ) $ 5 

3. . 10. .25.0.1.1.0 
0.0, .60,-1. ,1. 

CPCASE 

2 . 0. 3 . 0 

$MOST ADVERSE $ 

CPCASE 

3 . 0 . 5 . 0 
$FLAT T0P$ 

CPCASE 

4. 0. 2.0 


- a grid will be used 

- all three lines from same data file 

- first set of dependent data located in 
second coloumn of "CPCASE", use symbol 
3, no connecting lines 

- third column, symbol 5 


- fourth column, symbol 2 


Note: Data in this file must be in E20.5 format. 
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0.750 ALT = 38000 FT. 



Example 4 - continued 



Example 5 


This example illustrates the use of a new font using embedded text 
commands and the use of greater than one frame in a run. Also, the 
capability of the Metafile translator is shown on the third page. 


2 

$$ 

$$ 

$X/C$ 

$— C[BSUB] [BLC]P[ESUB] [ELC] 

$FLIGHT CONDITIONS $ 

$ (F3 . 1) $ 3 

$ (F3 . 1) $ 3 

2. . 2. .5.1. 1.1.1 
CPDATA1 

2. 0. 0.1 

$M1, C[BSUB] [BLC]L[ESUB] [ELC] 3$ 

CPDATA2 

2 , 0 , 0, 2 

$ M3, C [ BSUB ] [ BLC ] L [ ESUB ] [ ELC ] 4 $ 

$$ 

$$ 

$x/c$ 

$[ FONT=9 ] D [ FONT ]C[ BLC] [BSUB]P[ELC] [ESUB]$ - greek (font 9) delta (Ref. 2 
$ FLIGHT CONDITIONS $ 

$ (F3 . 1) $ 3 

$ (F4 . 2 ) $ 4 

1. . 2. .05.1.1.1.1 
DFDATA 

2. 0. 0.1 

$M1 , C(BSUB] [BLC]L[ESUB] [ELC] 3$ 


Metafile commands to create a page with two plots are as follows: 


S 

S 

s 

s 

D 

D 

D 


MF 1 META 

V 1 (-1 1 -1 1) 

V 2 (-1 0 -.65 .65) 

V 3 (0 1 -.65 .65) 

P 1 MF 1 V 1 

P 2 MF 1 V 1 

P 1 MF 1 V 2 P 2 MF 1 V 3 


- Default viewport 

- Left Viewport 

- Right Viewport 

- Plot 1 on one page 

- Plot 2 on one page 

- Plots 1 and 2 one one page 
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FLIGHT CONDITIONS 



Example 5 ■ continued 






Example 6 


This example illustrates the use of a logarithmic axis on a plot. 


1 

$$ 

$$ 

$S$ 

$Q$ 

$$ 

$ (F4 . 0) $ 4 

$ (F3 . 1) $ 3 

2 , 10 • , 0 • , 1 , 1 , 2,0 
0 . , 110 . ,. 1 , 1 . 
CSTAH 

2. 1.3.0 
CONEH 

2 . 1 . 4 . 0 


- vertical axis will be logarithmic 

- maxima and minima must be specified, 
vertical axis goes through a single 
cycle only 
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Example 7 


This example illustrates multiple frames, the use of a grid, symbols, 
titles and a legend. Also, subscripting is shown. 


3 

$$ 

$M = 0.750 ALT = 38000 FT. FREQ. 
$X/C$ 

$N [ BSUB ] CF [ ESUB ] $ 

$$ 

$ (F4 . 2 ) $ 4 

$ (F4 . 1) $ 4 

3. . 10.5.0.1.1.0 
0.0, .60,0.0,25. 

SALP016 

2. 0. 3.0 

$M0ST ADVERSE $ 

SALP017 

2. 0. 5.0 
$FLAT T0P$ 

SALP018 

2 . 0 . 2.0 

$LESS ADVERSE $ 

$$ 

$M = 0.750 ALT = 38000 FT. FREQ. 
$X/C$ 

$RE [ BSUB ] CF [ ESUB ] $ 

$$ 

$ (F4 . 2 ) $ 4 

$ (F6 . 1) $ 6 

3. . 10.50. .0.1. 1.0 
-.10, .60,-300. ,0.0 
MRCF016 

2. 0. 3.0 

$M0ST ADVERSE $ 

MRCF017 

2. 0. 5.0 
$FLAT T0P$ 

MRCF018 

2 . 0 . 2.0 

$LESS ADVERSE $ 

$$ 

$M = 0.750 ALT = 38000 FT. FREQ. 
$X/C$ 

$N[BSUB]TS [ESUB] $ 

$$ 

$ (F4 . 2 ) $ 4 

$ (F4 . 1) $ 4 

3, .10,5,0,1,1,0 
0.0,. 60 ,0.0,25. 

FILE1 


0.5 HZ. XLEN = 0.00065$ 


0.5 HZ. XLEN = 0.00065$ 


5000 HZ. XLEN = 0.007 
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Example 7 - continued 


2. 0. 3.0 

$MOST ADVERSE $ 
FILE2 

2. 0. 5.0 
$FLAT TOP$ 
FILE3 

2 . 0 . 2.0 

$LESS ADVERSE $ 
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ALT = 38000 FT. FREQ. = 0.5 HZ 










Example 8 


This example illustrates the use of symbols and lines in the same frame. 


1 

$FLIGHT DATA VS. THEORY $ 

$$ 

$X/C$ 

$-C [BSUB] [BLC] P[ESUB] [ELC] 

$$ 

$ (F3 . 1) $ 3 

$ (F3 . 1) $ 3 

4, .2, *5, 1,1, 1,0 
— .2, 1.0, -.5,2. 

JETCP4 

2 , 0 , 0,1 

$ THEORY-NO STEP$ 

D988 
2 , 0 , 2,0 
$FLT . 988$ 

RAMPCP8 

2.0. 0.3 
$THEORY-STEP$ 

D1112A 

2. 0. 3.0 
$FLT . 1070$ 


linestyle 1 
symbol 2 
linestyle 3 
symbol 3 
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Example 8 - continued 



Example 9 


This example shows the use of a line with selected symbols at each data 
point. A grid is also incorporated. An example of the data file is also 
shown . 


1 

$ SALES CHART $ 

$ [ F0NT=4 ] XTRA [ FONT ] MANUFACTURING COMPANY $ 
$MONTH$ 

$ NUMBER OF WIDGETS $ 

$YEAR$ 

$ (F3 . 0) $ 3 

$ (F5 . 1) $ 5 

3/1/20,0,1,1,1 


WIDDAT1 

2. 0. 5.1 
$1979$ 
WIDDAT1 

3. 0. 2.1 
$1981$ 
WIDDAT1 


4, 0,4,1 
$1983$ 


italics font; note: [FONT] 

returns to the default 
font 


use marker 5 with 
connecting lines 


The input file "WIDDATl” should have 
read from columns other than l and 2 

E20.5 format since 
. Example: 

data is to be 

. 10000E+01 

. 13500E+03 

. 10100E+03 

. 55000E+02 

. 20000E+01 

. 11200E+03 

. 87000E+02 

. 67000E+02 

. 30000E+01 

. 90000E+02 

. 76000E+02 

. 90000E+02 

. 40000E+01 

. 76000E+02 

. 55000E+02 

. 10400E+03 

. 50000E+01 

. 30000E+02 

. 90000E+02 

. 17500E+03 

. 60000E+01 

. 12000E+02 

. 91000E+02 

. 17800E+03 

. 70000E+01 

. 30000E+02 

. 14400E+03 

. 19000E+03 

. 80000E+01 

. 50000E+02 

. 16000E+03 

. 16000E+03 

. 90000E+01 

. 60000E+02 

. 16500E+03 

. 14300E+03 

. 10000E+02 

. 67 000E+02 

. 11200E+03 

. 90000E+02 

. 11000E+02 

. 50000E+02 

. 68000E+02 

. 35000E+02 

. 12000E+02 

. 13000E+02 

. 22000E+02 

. 34000E+02 
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MONTH 

Example 9 - continued 




Example 10 


This example shows an unusual line. 


1 

$C141 WING GEOMETRY $ 
$ETA = 0.683$ 

$X/C$ 

$Z/C$ 

$$ 

$ (F4 . 2 ) $ 4 

$ (F4 . 2 ) $ 4 

1. . 2. .25.0.1.1.0 
0. ,1. , - . 5 , .5 
C141 

2 . 0 . 0.1 

$ CURVE FIT$ 
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Appendix A 

Source Program Listing 
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ooo 


The following is a listing of the line plotting program source code. 
This is found in the file "LPLTR" under UN=826035N. 


PROGRAM LPLTR ( INPUT , OUTPUT , TAPE5=INPUT , TAPE8=0UTPUT ) 

CHARACTER*? FN, DFIL, PTITLE 
CHARACTER* 2 ISCPT,IMET 

CHARACTER* 5 0 LABE LX , LABELY , LEGEND , LEGN 
CHARACTER* 8 0 TITLE1 , TITLE2 
CHARACTER* 20 HFORM , VFORM 

INTEGER V , VS , LANS , LSTY , NFRAME , NLINES , MMAX , NPOINTS 
INTEGER TTYPE , IHLOG , I VLOG , HFW , VFW , FLAG 
DIMENSION V (3000) , DEPEND (7, 500) ,DEP(500) ,NPOINTS(7) ,IND(500) 

C 

C **NOTE: 

C DIMENSION V(VS) , DEPEND (LMAX, MMAX) ,DEP(MMAX) , NPOINTS (LMAX) ,IND(MMAX) 

C DIMENSION INDEP (LMAX, MMAX) , LEGN ( LMAX) , MARK ( LMAX ) , LAS (LMAX) 

C 

DIMENSION INDEP (7,5 00) , LEGN ( 7 ) ,MARK(7) , LAS (7) , PTITLE (7) 

REAL INDEP , DEPEND , HTICINC , VTICINC , IND , DEP 
INTEGER HAXIS , VAXIS , HTIC, VTIC 
INTEGER STR1 , STR2 , STR3 , TXIDX1 , TXIDX2 , QUAL 
COMMON /BLOCK/ NLINES 
DATA VS, MMAX, LMAX /3000,500,7/ 

DATA FLAG,NCHAR /0,1/ 

DATA STR1 , STR2 , STR3 , TXIDX1 , TXIDX2 , QUAL / 1 , 3, 5, 1,2, 4/ 

DATA IBOX9, SCX, SCY, TTYPE, LI /9 , 1 . , 1 . , 1 , 0/ 

PTITLE (1) = 'ONE' 

PTITLE ( 2 ) = 'TWO' 

PTITLE ( 3 ) = ' THREE ' 

PTITLE (4) = ! FOUR ' 

PTITLE (5) = 'FIVE' 

PTITLE (6) = 'SIX' 

PTITLE (7) = 'SEVEN' 

OBTAIN OUTPUT DESTINATION & CONTROL FILENAME FROM TERMINAL 

WRITE (8,*) ' PLOT TO SCREEN? (Y/N) ' 

READ (5,2222) ISCPT 
2222 FORMAT (A) 

WRITE (8,*) ' DO YOU WISH TO CREATE A METAFILE ? (Y/N) ' 

READ (5,2222) IMET 

WRITE (8,*) • ENTER CONTROL FILENAME: ' 

READ (5,10) FN 
10 FORMAT (A7) 

C 

C BEGIN DI3000 , GRAPHMAKER GRAPHICS 
C 

C INITIALIZATION - SET VIEWING WINDOW, LINESTYLES, TEXT 
C 

IF (ISCPT .EQ. *Y '.OR. ISCPT .EQ.' Y') THEN 
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668 


60 


669 

C 

C 

C 


41 
40 

42 

43 


C 

C 

C 


C 

C 


CALL JCHINI ( . TRUE . , 1 ) 

ELSE 

CALL JCHINI ( . TRUE . , 0 ) 

IMET = 1 Y ' 

GO TO 668 
END IF 

IF (IMET .EQ. 'N ' .OR. IMET .EQ. ' N') GO TO 668 
CALL JDINIT ( 0) 

CALL JDEVON ( 0 ) 

CALL JVSPAC(-1. ,1. ,-.7, .7) 

CALL J CHART (V , VS) 

CALL JCHEXT(V,VS, 0,1000. ,700. ) 

CALL JXTEXT (V, VS , TXIDX1 ,3, 0.0,1. 4,0) 

CALL JXTEXT (V, VS , TXIDX2 ,3, 0.0, 1.0,0) 

CALL JXBOX (V, VS , IBOX9 ,8,0, .FALSE.) 

LSTY=128 

DO 60 III =1,9 

CALL JXLINE (V, VS , III , 0 , 16383 , LSTY, 20000) 

LSTY = LSTY + 1 
CONTINUE 

CALL JXLINE (V, VS, 10, 0,16383, 0,28000) 

CALL JXLINE (V, VS, 11, 0,16383, 0,10000) 

OPEN (UNIT=4,FILE=FN) 

IF (IMET .EQ. 'N ' .OR. IMET .EQ. ' N') GO TO 669 
CALL JMSTRG ( 0 , PTITLE ( 1 ) ) 

CONTINUE 

READ CONTROL DATA 

READ (4,*) NFRAME 
DO 200 I = 1, NFRAME 
READ (4,42) TITLE1 
READ (4, 42) TITLE2 
READ (4,43) LABE LX 
READ (4,43) LABELY 
READ (4,43) LEGEND 
READ (4, 41) HFORM , HFW 
READ (4,41) VFORM , VFW 
FORMAT (A10, 14) 

FORMAT (A30) 

FORMAT (A80) 

FORMAT (A50) 

READ (4,*) NLINES , HTICINC , VTICINC , IGRID , IHLOG , IVLOG , I INP 
IF (NLINES . GT . LMAX) THEN 

WRITE ( 8 , * ) '***> INCREASE LMAX AND DIMENSIONS !' 

STOP 
END IF 

IF IINP=0 THEN READ INPUT FOR MAXIMA & MINIMA 

IF (IINP.EQ.O) THEN 

READ (4,*) HMIN , HMAX , VMIN , VMAX 
FLAG = 1 

DEFAULT TIC INCREMENTS 
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c 

IF (HTICINC .EQ. 0.) THEN 
HTICINC= (HMAX-HMIN) /10 
END IF 

IF ( VTICINC . EQ . 0 . ) THEN 
VTICINC= (VMAX-VMIN) /10 
END IF 
ENDIF 

DO 150 J = 1 , NLINES 
READ (4,10) DFIL 
OPEN (UNIT=6 , FILE=DFIL) 

READ ( 4 , * ) NCOL, LANS, MARK ( J) , LAS (J) 

IF (LANS.NE. 0) GOTO 2 
C 

C LEGEND ENTRY 
C 

READ (4, 43) LEGN(J) 

C 

C READ DATA FROM PROPER COLUMN 
C 

2 IF (NCOL.EQ.2) THEN 

DO 11 M=1 , MMAX 

READ ( 6 , * , END=2 3 ) INDEP ( J , M) , DEPEND ( J , M) 

12 FORMAT (2E15. 5) 

11 CONTINUE 

ENDIF 

IF (NCOL. EQ. 3 ) THEN 
DO 14 M=1 , MMAX 

READ (6,15, END=23 ) INDEP (J, M) , DEPEND (J,M) 

15 FORMAT ( E2 0.5,20X,E20.5) 

14 CONTINUE 

ENDIF 

IF (NCOL.EQ.4) THEN 
DO 17 ,M=1 ,MMAX 

READ (6,18, END=23 ) INDEP ( J, M) , DEFEND (J,M) 

18 FORMAT (E20 . 5 , 40X, E20 . 5) 

17 CONTINUE 

ENDIF 

IF (NCOL. EQ. 5) THEN 
DO 20 M=1 , MMAX 

READ (6, 21,END=23) INDEP (J ,M) , DEPEND (J , M) 

21 FORMAT (E20 . 5, 60X, E20 . 5) 

20 CONTINUE 

ENDIF 

IF (NCOL. EQ. 6) THEN 
DO 25 M = 1,MMAX 

READ (6,24, END=2 3 ) INDEP(J,M) , DEPEND ( J,M) 

24 FORMAT (E20.5,80X,E20.5) 

25 CONTINUE 
END IF 

IF (NCOL.GT. 6) THEN 

WRITE (8,*) 'YOU WILL HAVE TO MODIFY THE PROGRAM! ! ! ! ' 
STOP 
ENDIF 
C 
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c 

c 

23 

150 

C 

C 

C 

C 


667 

C 

C 

C 

C 


C 

C 

C 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


NUMBER OF POINTS IN DATA FILE IS NOW KNOWN 

NPOINTS(J) = M - 1 
REWIND (UNIT=6) 

CONTINUE 


IF MAXIMA AND MINIMA HAVE ALREADY BEEN SET, BEGIN PLOT 
OTHERWISE FIND THEM USING SUBROUTINE MAXFNDR & DATA 


IF ( FLAG . EQ . 1 ) GO TO 667 

CALL MAXFNDR ( INDEP , NPOINTS , HTICINC , HMAX , HMIN) 
CALL MAXFNDR (DEPEND, NPOINTS , VTICINC , VMAX , VMIN) 
IF (I .EQ. 1) CALL JGRAPH (V, VS , I) 


800 X 600 IS THE SIZE OF THE DATA SPACE IN CHART UNITS 
40, 150 IS THE POSITION OF THE LOWER LEFT HAND CORNER 


CALL JDSWDW (V, VS , NCHAR, 800., 600.) 
CALL JDSPOS (V, VS , NCHAR, 40., 150.) 
CALL JTXQAL ( V , VS , QUAL) 

CALL JTXHGT(V,VS, 35. 0,1. 0,1.0) 
CALL JTXBOX (V, VS ,0,0, TXIDX1) 

TITLE1 


CALL JSTNOT ( V , VS , NCHAR, STR1 , TITLE 1 ) 
CALL JPONOT (V, VS, NCHAR, STR1, 500. ,890. ) 
CALL JTXBOX (V, VS, 0,0, TXIDX1) 

TITLE2 


CALL JSTNOT (V, VS, NCHAR, STR3,TITLE2) 
CALL JPONOT(V, VS, NCHAR, STR3, 500. ,820. ) 
CALL JTXBOX (V, VS, 0,0, TXIDX1) 

LEGEND 


CALL JTXHGT (V, VS , 30 . , 1 . , 1 . ) 

CALL JSTNOT ( V , VS , NCHAR , STR2 , LEGEND) 
CALL JPONOT (V, VS, NCHAR, STR2, 990. ,720. ) 
CALL JTXBOX (V, VS, 0,0,TXIDX2) 

HORIZONTAL AXIS 


HAXIS=1 

CALL JTXHGT (V, VS, 30. 0,1. 0,1.0) 

CALL JSTHAX ( V , VS , NCHAR, HAXIS , HMIN , HMAX , LABE LX) 
CALL JAXATR ( V , VS , NCHAR, HAXIS , 10 , SCX , IHLOG) 

VERTICAL AXIS 


VAXIS=2 

CALL JSTVAX (V, VS , NCHAR, VAXIS , VMIN , VMAX , LABELY) 
CALL J AXATR ( V , VS , NCHAR , VAXIS , 1 0 , S C Y , IVLOG ) 
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C TICK MARKS/GRID LINES & LABELS 
C 

C HORIZONTAL TICKS/GRID 

CALL JTXBOX (V, VS ,0,0, TXIDX2 ) 

CALL JTXHGT(V,VS, 25. 0,1. 0,1.0) 

HTIC=1 

IF (IGRID.EQ. 0) TTYPE=3 

CALL JTIC ( V , VS , NCHAR , HAXIS , HTIC , HMIN , HMAX , HTI CINC ) 

CALL JTCTYP (V, VS , NCHAR, HAXIS , HTIC , TTYPE , 10) 

IF ( IGRID . EQ . 1 ) CALL JTCATR (V, VS, NCHAR, HAXIS, HTIC, -10. , 0. , 10) 
CALL JSTFMT (V, VS , NCHAR, HAXIS , HTIC , HFW, HFORM) 

C 

C MINOR TICKS 
C 

HTIC=2 

HTICINC=HTICINC/5 
TTYPE = 2 

CALL JTIC ( V , VS , NCHAR , HAXIS , HTIC , HMIN , HMAX , HTICINC ) 

CALL JTCTYP(V, VS, NCHAR, HAXIS, HTIC, TTYPE, 11) 

CALL JTCATR (V, VS, NCHAR, HAXIS, HTIC, -5. ,0. ,11) 

C 

C VERTICAL TICKS/GRID 
C 

300 CALL JTXJST(V,VS, 3,2) 

TTYPE=1 

VTIC=1 

CALL JTIC ( V , VS , NCHAR , VAXIS , VTIC , VMIN , VMAX , VTICINC ) 

IF (IGRID.EQ. 0) TTYPE=3 

CALL JTCTYP ( V , VS , NCHAR , VAXIS , VTIC , TTYPE , 10 ) 

IF ( IGRID . EQ . 1 ) CALL JTCATR(V, VS , NCHAR, VAXIS , VTIC , -10 ., 0 ., 10) 
CALL JSTFMT ( V , VS , NCHAR , VAXIS , VTIC , VFW , VFORM) 

C 

C MINOR TICKS 
C 

VTIC=2 

VTICINC=VTICINC/5 

TTYPE=2 

CALL JTIC ( V , VS , NCHAR , VAXIS , VTIC , VMIN , VMAX , VTICINC ) 

CALL JTCTYP (V, VS, NCHAR, VAXIS, VTIC, TTYPE, 11) 

CALL JTCATR (V, VS, NCHAR, VAXIS, VTIC, -5. ,0. ,11) 

301 CONTINUE 
C 

C PASS DATA 
C 

DO 75 II = 1 , NLINES 

DO 72 JJ = 1 , NPOINTS (II) 

IND ( JJ) =INDEP (II , JJ) 

DEP(JJ) = DEPEND ( II, JJ) 

72 CONTINUE 

IIP1=(II-1) *2+1 

CALL JRDATA (V, VS , IIP1, IND, NPOINTS (II) ) 

CALL JINDEP (V, VS, NCHAR, IIP1) 

IIP1=IIP1+1 

CALL JRDATA (V, VS, IIP1 , DEP, NPOINTS (II) ) 

CALL JDEPEN (V, VS, NCHAR, II, IIP1) 
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OUO UOU OUU HUUU 


LINES OR SYMBOLS 

IF(MARK(II) . EQ. 0 . AND. LAS (II) .EQ.O) LAS(II) = 1 
CALL JDTATR(V, VS , NCHAR, II , 0 , MARK (II) , LAS (II) ) 
LI = LISTOR 


LEGEND 

IF (LANS .NE . 0) GO TO 101 
CALL JTXBOX (V, VS , 0 , 0 , TXIDX1) 

CALL JSTLGD (V, VS, NCHAR, ’$$' ) 

POSITION OF THE LEGEND 

PUT = 680. - II * 40. 

CALL JLGPOS (V, VS, NCHAR, 870. , PUT) 

CALL JTXHGT (V, VS ,27.0,1.0,1.0) 

CALL JSDLGD (V, VS , NCHAR, II , LEGN (II) ) 

1 CONTINUE 

SHOW CHART 

CALL JCHSHW (V, VS , — 1 . ,1. ,-.7, .7) 

IF (II .LT. NLINES) GO TO 201 

IF ( ISCPT .EQ. 'Y ' .OR. ISCPT .EQ. ' Y') CALL JPAUSE(l) 
IF (I .GE. NFRAME) GO TO 200 
IF (IMET .EQ. 'N ' .OR. IMET .EQ. ' N') GO TO 199 
CALL JMSTRG (0 , PTITLE (1+1) ) 

199 CALL J FRAME 

201 CONTINUE 

CALL J CHART (V, VS) 

CALL JGRAPH(V, VS, NCHAR) 

CALL JCHEXT (V, VS, 0,1000. ,700. ) 

CALL JXTEXT (V, VS , TXIDX1 ,3, 0.0, 1.4,0) 

CALL JXTEXT (V, VS ,TXIDX2 ,3, 0.0, 1.0,0) 

CALL JXBOX ( V, VS , IBOX9 ,8,0,. FALSE . ) 

IF (II .LT. NLINES) THEN 

CALL JDSWDW (V, VS , NCHAR, 800., 600.) 

CALL JDSPOS (V, VS , NCHAR, 40 . , 150 . ) 

CALL JTXQAL ( V, VS , QUAL) 

CALL JSTHAX(V, VS, NCHAR, HAXIS,HMIN,HMAX, *$$') 

CALL JSTVAX(V, VS, NCHAR, VAXIS,VMIN, VMAX, '$$') 

CALL JAXATR ( V , VS , NCHAR , HAXI S , 1 0 , SCX , IHLOG ) 

CALL JAXATR (V, VS, NCHAR, VAXIS, 10, SCY,IVLOG) 

HM1 = HMIN - 10. 

VM1 = VMIN - 10. 

CALL JTIC ( V , VS , NCHAR , VAXIS , VTIC , VM1 , VM1 , 1 . ) 

CALL JTIC ( V , VS , NCHAR , HAXIS , HTIC , HM1 , HM1 , 1 . ) 

END IF 

LSTY=128 

DO 61 III =1,9 

CALL JXLINE(V, VS, III, 0, 16383, LSTY, 20000) 

LSTY = LSTY + 1 
61 CONTINUE 
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CALL JXLINE(V,VS, 10,0, 16383,0,28000) 

CALL JXLINE(V,VS, 11, 0,16383 ,0,10000) 

75 CONTINUE 

200 CONTINUE 

CALL JCHTRM ( . TRUE . ) 

STOP 

END 

C 

C THIS SUBROUTINE FINDS THE MINIMA AND MAXIMA OF A 2-D DATA 
C SET GIVEN THE NUMBER OF POINTS, AND THE DESIRED INCREMENT 
C BETWEEN DATA TIC MARKS SUCH THAT THE CLOSEST EVEN MAX AND 
C MIN ARE FOUND. 

C 

SUBROUTINE MAXFNDR (ARRAY , NPOINTS , INC , MAX , MIN) 

DIMENSION ARRAY (7,500) , NPOINTS ( 7 ) 

REAL MAX, MIN, INC 
COMMON /BLOCK/ NLINES 
MIN=ARRAY (1,1) 

MAX= ARRAY (1,1) 

598 DO 600 II = 1, NLINES 

DO 550 JJ = 1, NPOINTS (II) 

IF (MIN .GT. ARRAY (II , JJ) ) MIN=ARRAY ( I I , J J ) 

IF (MAX .LT. ARRAY (II , JJ) ) MAX = ARRAY(II,JJ) 

550 CONTINUE 

600 CONTINUE 

IF (MAX.GE.0) THEN 
DO 10 IJ = 1, 100 
TEMP=INC*IJ 

IF (TEMP .GE. MAX) THEN 
DIFF = ABS (TEMP-MAX) 

IF (DIFF. LE. INC) THEN 
MAX = TEMP 
GO TO 20 
ELSE 

WRITE (8,*) * NO EVEN MAX WAS FOUND' 

GO TO 20 
END IF 
END IF 

10 CONTINUE 
ELSE 

DO 15 IJ = 1, 100 
TEMP=— INC*IJ 
IF (TEMP .LT. MAX) THEN 
TEMP1=- ( I J - 1 ) *INC 
DIFF = ABS (TEMP1-MAX) 

IF (DIFF. LE. INC) THEN 
MAX = TEMPI 
GO TO 20 
ELSE 

WRITE (8,*) 'NO EVEN MAX WAS FOUND' 

GO TO 20 
END IF 
END IF 

15 CONTINUE 
END IF 
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20 CONTINUE 

TEMP=MAX 
DO 30 IJ = 1,100 
TEMP=TEMP- INC 
IF (TEMP. LE. MIN) THEN 
DIFF = ABS (MIN-TEMP) 

IF (DIFF. LE. INC) THEN 
MIN=TEMP 
GO TO 40 
ELSE 

WRITE (8,*) 'NO EVEN MIN WAS FOUND' 
GO TO 40 
END IF 
END IF 
30 CONTINUE 

40 CONTINUE 

RETURN 
END 
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Appendix B 
Notes to the User 
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Notes to the User 


This line plotting program was developed before the release of the 
Common Graphics Library (CGL) . Hence, it should be noted that a 
modification of this program might be to include calls to CGL 
subroutines in order to obtain the standard NASA markers and 
linestyles on plots. Color was not implemented in order to keep 
the program device independent. However, by using the color option 
(which is presently set to zero) in calls to subroutine JXLINE, 
lines can be distinguished using color on terminals which have that 
capability. 

Also, to simplify the process of obtaining plots, the length of the 
axes is fixed for all plots. To scale axes using this program, 
user input axes lengths need to be added to the control file 
format. Then calls to JVSPAC, JCHEXT, JDSWSW, JDSPOS, and JPONOT 
must be modified to set the view space, the chart extent, the data 
space, the position of the lower left hand corner of the data space 
and the positions of the plot text, respectively, in terms of the 
input axes lengths. JVSPAC requires the ratio of the shorter axis 
to the longer axis. JCHEXT defines the picture coordinate system 
and all other calls requiring coordinate information refer to this 
system. Whether the chart extent information (JCHEXT) is input 
directly from the control file or calculated from data space inputs 
(JDSWDW) is up to the user, but parameters for all the other 
subroutines may be calculated from that information (Ref. 1) . 

Finally, since data files containing more than one set of dependent 
variables are read in E20.5 format, it might be useful to customize 
the program by either changing the FORMAT statements in lines 12, 
15, 18, 21 and 24, or by adding an input to read the data format 
from the control file. 
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